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Abstract. BL Lacertae has been the target of four observing campaigns by the Whole Earth Blazar Telescope (WEBT) 
collaboration. In this paper we present UBVRI light curves obtained by the WEBT from 1994 to 2002, including the 
last, extended BL Lac 2001 campaign. A total of about 7500 optical observations performed by 31 telescopes from 
Japan to Mexico have been collected, to be added to the ~ 15600 observations of the BL Lac Campaign 2000. All 
these data allow one to follow the source optical emission behaviour with unprecedented detail. The analysis of the 
colour indices reveals that the flux variability can be interpreted in terms of two components: longer-term variations 
occurring on a few-day time scale appear as mildly-chromatic events, while a strong bluer-when-brighter chromatism 
characterizes very fast (intraday) flares. By decoupling the two components, we quantify the degree of chromatism 
inferring that longer-term flux changes imply moving along a ~ 0.1 bluer- when-brighter slope in the B — R versus R 
plane; a steeper slope of ~ 0.4 would distinguish the shorter-term variations. This means that, when considering the 
long-term trend, the i?-band flux level is related to the i?-band one according to a power law of index ~ 1.1. Doppler 
factor variations on a "convex" spectrum could be the mechanism accounting for both the long-term variations and 
their slight chromatism. 



1. Introduction 

Blazars are a class of radio-loud active galactic nuclei 
(AGNs) well known for their peculiar properties, such as 
their intense and extremely variable non-thermal emis- 
sion across all the electromagnetic spectrum from the 
radio band to 7-rays (sometimes up to TeV energies). 
Blazars are divided into two subclasses: flat-spectrum ra- 
dio quasars and BL Lac objects, whose main difference is 
the lack or weakness of emission lines in BL Lac objects. 

The commonly accepted paradigm capable of explain- 
ing most of the observed radio-loud AGN properties in- 
volves a supermassive black hole surrounded by an accre- 
tion disc, feeding a relativistic plasma jet emitting syn- 
chrotron and inverse-Compton radiation. In this scenario 
blazars would be the fraction of sources with the jet ori- 
ented at a small angle with respect to the line of sight. 

To investigate the details of the physical processes and 
geometric conditions at the base of the extreme emis- 
sion variability, a great observing effort is needed, pos- 
sibly organizing multiwavelength campaigns to follow the 
source emission behaviour simultaneously at different fre- 
quencies. For example, the Whole Earth Blazar Telescope 
(WEBT; http://www.to.astro.it/blazars/webt/, e.g. 
Villata et al. 12170771 127701)1 is an international collabora- 
tion of optical and radio observers whose aim is the inten- 
sive and accurate monitoring of selected blazars during 
time-limited campaigns. The distribution in longitude of 
the WEBT members makes continuous (24 hours per day) 
optical monitoring possible, at least in principle. Since its 
inception in 1997, various BL Lac objects have been the 
target of the WEBT 1 ; in particular, four campaigns have 
been devoted to BL Lacertae, the prototype of the above 
mentioned subclass. 

The first two BL Lac campaigns were organized in 
1999 in conjunction with pointings of the X-ray satel- 
lites BeppoSAX (June 5-7) and ASCA (June 28-30). 
They lasted a few days only and involved a restricted 
participation of the WEBT members. The results of 
the BeppoSAX- WEBT campaign have been published in 
Ravasio et al. (2003): the shape of the X-ray spectrum 

1 Targets of WEBT campaigns have been, besides BL 
Lacertae, AO 0235+16 (Raiteri et al. HTj7TJ, Mkn 421, S5 
0716+71 (Villata et aL l^DODll . Mkn 501, 3C 66A. 



was concave, with a very hard component above 5-6 keV; 
moreover, a very fast variability event was detected at the 
lowest X-ray energies, with no simultaneous optical coun- 
terpart. 

The subsequent WEBT BL Lac 2000 campaign was an 
impressive observing effort with a long extension in time: 
more than 15000 optical observations were performed by 
24 telescopes in 11 countries from May 2000 to January 
2001. The results of these observations are reported in 
Villata et al. (2002J). The original motivation for the cam- 
paign was to provide the low-energy observing counter- 
part for a high-energy campaign organized in July-August 
2000, involving X-ray satellites as well as TeV detectors 
(Bottcher et al. I2003fl . In this period (the "core cam- 
paign") the densest sampling was achieved, with observ- 
ing gaps limited to a few hours, mainly due to the lack of 
observers in the Pacific area. 

After the BL Lac Campaign 2000, the WEBT collab- 
oration launched a new BL Lac campaign in June 2001 
(due to the detection of a fast brightening of the source), 
which lasted until the end of February 2002. In this paper 
we present the optical light curves of this BL Lac 2001 
campaign and of its so-called "extension" , which includes 
the two WEBT campaigns of 1999 and a collection of un- 
published data back to 1994. We also include in our anal- 
ysis previously published optical data from the BL Lac 
Campaign 2000 (Villata et al. 12002(1 and literature data 
from 1994 to 2002. As a result, the optical light curves 
that we present in this work, based on new and previously 
published WEBT and literature data, cover a period of 
more than 8 years (9 observing seasons). Radio and opti- 
cal light curves covering the period 1968-2003 will be pre- 
sented in a forthcoming paper (Villata et al. 2004), where 
timing and cross-correlation analyses will be performed on 
both optical and radio light curves. 

In Sect. 2 we report the optical observations of BL 
Lacertae carried out by the WEBT collaboration; the 
UBVRI light curves from 1994 to 2002 are presented in 
Sect. 3, and a detailed colour analysis is performed in Sect. 
4. Conclusions are drawn in Sect. 5. 
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Table 1. List of participating observatories by longitude: d is the telescope diameter; N Q ^ S is the total number of 
observations done, i.e. the number of unbinned data; Nu, Nb, Ny, Nr, and Nj are the numbers of data points in 
UBVRI remained after discarding and binning some of the original data; the numbers in brackets refer to the WEBT 
campaign 2001 alone (May 2001 - February 2002). 



Observatory 


d (cm) 


N ohs 


N v 


N B 


N v 


N R 


Nj 


Kyoto, Japan 


25 


190 (0) 











19 (0) 





Yunnan, China 


100 


25 (25) 





2 (2) 


2 (2) 


7 (7) 





Mt. Maidanak (T-60), Uzbekistan 


60 


50 (50) 


15 (15) 


16 (16) 


16 (16) 








Mt. Maidanak (AZT-14), Uzbekistan 


50 


150 (150) 


49 (49) 


47 (47) 


49 (49) 


3 (3) 





Mt. Maidanak (AZT-22), Uzbekistan 


150 


418 (418) 


50 (50) 


87 (87) 


65 (65) 


108 (108) 


61 (61) 


Abastumani, Georgia (FSU) 


70 


1823 (316) 





276 (36) 


232 (0) 


903 (213) 


150 (0) 


Crimean, Ukraine 


70 


267 (267) 





67 (67) 


67 (67) 


67 (67) 


66 (66) 


Skinakas, Crete 


130 


510 (239) 





186 (119) 


69 (0) 


184 (118) 


67 (0) 


Tuorla, Finland 


103 


80 (0) 








79 (0) 








MonteBoo, Czech Republic 


62 


754 (754) 











282 (282) 





Vallinfreda, Italy 


50 


98 (0) 








49 (0) 








Perugia, Italy 


40 


524 (419) 








112 (112) 


208 (180) 


166 (125) 


Greve, Italy 


32 


55 (0) 














55 (0) 


Torino, Italy 


105 


487 (54) 





104 (1) 


70 (2) 


271 (24) 


12 (2) 


Sabadell, Spain 


50 


32 (32) 








3 (3) 


12 (12) 





Guadarrama, Spain 


20 


3 (3) 








1 (1) 


1 (1) 





Roque (KVA), La Palma 


60 


322 (0) 











152 (0) 





Roque (KVA), La Palma 


35 


42 (42) 








10 (10) 


10 (10) 





Roque (JKT), La Palma 


100 


64 (64) 





9 (9) 


10 (10) 


32 (32) 


6 (6) 


Roque (NOT), La Palma 


256 


44 (44) 


4(4) 


7(7) 


7(7) 


7(7) 


3(3) 


Roque (TNG), La Palma 


358 


138 (0) 


3(0) 


62 (0) 


5(0) 


62 (0) 


2(0) 


Olin, Connecticut 


51 


400 (361) 








8(0) 


78 (64) 


85 (71) 


Hopkins, Massachusetts 


60 


12 (0) 














11 (0) 


Foggy Bottom, New York 


40 


248 (83) 











248 (83) 





Boltwood, Canada 


40 


96 (96) 





6(6) 


6(6) 


6(6) 


5(5) 


Bell, Kentucky 


60 


23 (23) 











22 (22) 





St. Louis, Missouri 


36 


423 (340) 





10 (4) 


40 (16) 


161 (137) 


39 (15) 


Sommers-Bausch, Colorado 


60 


9(9) 








4(4) 


4(4) 





Lowell, Arizona 


180 


36 (0) 








18 (0) 





18 (0) 


San Pedro Martir, Mexico 


210 


9(9) 





3(3) 


3(3) 


3(3) 





San Pedro Martir, Mexico 


150 


123 (123) 


4(4) 


12 (12) 


11 (11) 


11 (11) 


4(4) 


Total 




7455 (3921) 


125 (122) 


894 (416) 


936 (384) 


2861 (1394) 


750 (358) 



figurcl.jpg 



Fig. 1. UBVRI light curves of BL Lacertae from 1994 to 2002; the numbers on the left indicate the number of data 
points in the extended BL Lac Campaign 2001 (blue/dark), in the BL Lac Campaign 2000 (red/grey), and from 
the literature (light cyan/light grey); yellow (shaded) strips correspond to the periods of the 1999 and 2001 WEBT 
campaigns. 



2. Observations of BL Lacertae by the WEBT 
collaboration 

The BL Lac Campaign 2001 was started in June 2001 
after the detection of a fast brightening of the source. 



Optical and radio monitoring was continued until the end 
of February 2002. However, at the time of data analysis, 
it was decided to collect all optical and radio data taken 
by the members of the WEBT collaboration dating back 
to 1968. The observations performed out of the proper 
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Fig. 2. Folded i?-band light curve of BL Lacertae from 1997 to 2002; symbols as in Fig.Q] 



campaign period (May 2001 - February 2002) constitute 
the so-called "campaign extension" ; they also include the 
still unpublished data of the WEBT campaign of late June 
1999 and the partially published (Ravasio et al. 2002) data 
of the early- June 1999 campaign. In this first paper, we 
present the optical data only, starting from 1994. 

Table 1 shows the list of optical observatories that par- 
ticipated in the data collection with unpublished data or 
data from the early- June 1999 campaign. The name and 
location of the observatory (Col. 1) is followed by the 
telescope diameter (Col. 2), by the total number of ob- 
servations made (Col. 3) and the number of data points 
in UBVRI derived, sometimes after binning (Cols. 4-8). 
The numbers in brackets refer to the WEBT campaign 
2001 alone. 

The observing strategy during the BL Lac Campaign 
2001 was nearly the same as in the BL Lac Campaign 2000 
(Villata et al. l2()()2|l : in general, observers with telescopes 
larger than 60 cm were invited to perform BR sequences 
of frames during each observing night, along with a com- 
plete (U)BVRI sequence at the beginning and at the end. 
Exposure times were chosen to obtain a good compro- 
mise between high precision and high temporal density. 
Participants with smaller-size telescopes were suggested 
to carry out observations in the R band only. 

In general, CCD data were collected as instrumental 
magnitudes of the source and reference stars to apply the 
same analysis and calibration procedures to all datasets. 
Instrumental magnitudes were obtained by either stan- 
dard aperture photometry procedures or Gaussian fitting, 
in most cases by the observers themselves. A minority of 
data came from photometer observations: these data were 
provided directly as standard magnitudes. 

A careful data assembly was required because of the 
huge amount of data coming from many different instru- 
ments: details can be found in Villata et al. ( 2002). 



3. Optical light curves 

UBVRI light curves from 1994 to 2002 are shown in Fig. 
^ Blue (dark) dots refer to data belonging to the BL Lac 
Campaign 2001 plus extension (see Sect. 2); red (grey) 
points represent the BL Lac Campaign 2000 data (Villata 



et al. 12002(1 : literature data 2 are shown as light cyan (light 
grey) dots. The periods of the 1999 and 2001 WEBT cam- 
paigns on BL Lacertae are highlighted by yellow (shaded) 
strips. 

Each panel in the figure covers a range of 4 magni- 
tudes. One can easily see that the variation amplitude 
is greater at higher frequencies, which is a common be- 
haviour in blazars. For instance, during the well-sampled 
1997 outburst, the brightness excursion was 3.09 mag in 
the B band, 3.06 mag in V, 2.94 mag in R, and 2.47 mag 
in /. Hence, the variation in the B band was 5% wider 
than in the R band and 25% wider than in the / band. 
Actually, these values of the AB excess have to be in- 
tended as underestimates, because the _B-band light curve 
is undersampled, in particular around maxima and min- 
ima. In general, in this period the number of B data is 
only 18% of the R one and 44% with respect to the I 
data. 

Another example can be derived by considering the 
period of the WEBT campaign 2001. Total magnitude 
variations were AU = 2.63, AB = 2.34, AV = 2.32, 
AR = 2.12, and Al = 2.07. Thus, again in spite of their 
poorer sampling (see Table 1), the UBV data show a 
larger total variation with respect to the R (and /) band. 

From Fig. n one can also notice that, while the mini- 
mum brightness states are very different, there seems to be 
a well-defined brightness upper limit. Indeed, the bright- 
est levels detected in summer 1997, late 1999, mid 2001 
and early 2002 reached about the same value. 

The best-sampled light curve is the i?-band one (8672 
data points); in Fig. [21 we show the most interesting pe- 
riod, 1997-2002, folded in three panels. Each panel covers 
two observing seasons, starting and ending on February 
17 (solar conjunction). The most noticeable features are 
the 1997 outburst, well covered by literature data, and 

2 Literature data have been taken from Maesano et al. 
I1997H . Fan et al. 119981 with dates taken from Fan & Lin 
HEOJl, Webb et al. (IT^Hl . Bai et al. itT^l . Tosti et al. 
03551 . Ghosh et al. JSPOot . Katajainen et al. (ISuTTOt . Clements 
& Carini (J5D5TJ|, Fan et al. (ISOOlTl . Sobrito et al. (|2770H . 
Papadakis et al. ( 2003 1 , Tosti et al. (12U04I . Sometimes the data 
used in this paper do not exactly correspond to the original 
ones since some correction/removal (often done in agreement 
with the data owners) was needed. Moreover, some datasets 
were completely re-analysed. 
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Fig. 3. Enlargements of the BL Lacertae i?-band light curve corresponding to the WEBT campaigns of late June - 
early July 1999 (top panel) and May 2001 - February 2002 (middle and bottom panels); symbols as in Fig. QJ in the 
top panel the yellow (grey) strip indicates the period of the ASCA pointing. 



the other two strong outbursts occurring in mid 2001 and 
early 2002, during the BL Lac 2001 WEBT campaign. 
Other well-sampled periods are mid 1999 (especially the 
WEBT campaigns) and of course the BL Lac Campaign 
2000. The 1997 outburst still remains the most spectac- 
ular event for both its magnitude excursion during the 



dimming phase and the large amplitude of the short-term 
variations (see e.g. Villata et al. 2002 for references). 

The three panels in Fig. [3] display the details of the 
WEBT campaigns carried out in mid 1999 contempora- 
neously with the ASCA pointing (top) and in 2001-2002 
(middle and bottom); they correspond to the second and 
third yellow (shaded) areas of Fig. [21 Details of the 1999 
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Fig. 4. An enlargement of the i?-band light curve of BL Lacertae on June 29, 1999, corresponding to the green (grey) 
box in Fig. the yellow (grey) filled circles represent the simultaneous i?-band light curve shifted by —1.62 mag. 



figure5.jpg 



Fig. 5. B — R (top panel) and V — I (middle panel) colour indices versus time from 1997 to 2002 compared to the 
i?-band light curve; the host-galaxy contribution has been removed from the data; symbols as in Fig.^ the horizontal 
lines divide each panel into four equal ranges (1 mag for the light curve and 0.25 mag for the colour indices) for clarity. 



WEBT-BeppoSAX campaign (first yellow/shaded strip 
in Fig. [3$ have already been presented in Ravasio et al. 
(2002). In the upper panel of Fig. [21 spanning the period 
June 27.5 - July 3.5, 1999, the time coverage by the ASCA 
satellite is indicated by a yellow (grey) strip; a comparison 
between the optical and X-ray data will be performed else- 
where. A very well-sampled intranight variation was ob- 
served at the beginning of June 29 (green/grey box), which 
is better displayed in Fig. 2J where both the i?-band and 
the B-band data are shown to reveal the expected larger 
variation in the B band: AB = 0.243 versus AR = 0.219 
in 2.5 h. Notice how data from different telescopes are in 
fair agreement and how the well-defined trend confirmed 
by high-precision data in two bands provides the finest 
details of the variation. 



In the 2001-2002 WEBT campaign, one can notice the 
double-peaked, nearly symmetric outburst of mid 2001, 
with fast oscillations superimposed (Fig. |3 middle panel) . 
On the contrary, the outburst occurring at the beginning 
of 2002 appears as a peculiar event with respect to the 
previous ones because of its single, short-duration peak. 
However, one can recognize fast (intranight) flares per- 
turbing the main trend (Fig. |3 bottom panel). 

4. Colour analysis 

As done for the BL Lac Campaign 2000 (Villata et al. 
2002$, the first step of the colour analysis is the subtrac- 
tion of the host-galaxy contribution from the fluxes to 
avoid its colour contamination. We have estimated the 
host-galaxy magnitudes in the various bands using data 
published by Scarpa et al. ((2000) and Mannucci et al. 
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Fig. 6. An enlargement of the B — R versus time curve (top panel) compared with the i?-band light curve (bottom 
panel) during the best-sampled period of the BL Lac Campaign 2001; orange (grey) dots represent the data used for 
calculating the colour indices. 



(2002), obtaining C7 host = 17.65, B host = 17.15, V host = 
16.16, -Rhost = 15.55, and /host = 14.92. Conversion of 
magnitudes into fluxes has been done as in Villata et al. 
( 2002J, where it was also estimated that only about 60% of 
the host-galaxy flux must be subtracted from the original 
data. 

Once this procedure was performed, in BL Lac 2000 
the resulting colour indices showed a fairly clear be- 
haviour: long-term variations appeared essentially achro- 
matic, while fast flares exhibited a strong bluer- when- 
brighter trend. 

In Fig. we show the time evolution of the B — R 
and V — I BL Lac colour indices from 1997 to 2002, com- 
pared with the galaxy-subtracted i?-band light curve. We 
have now a dataset which is much more extended both in 
time and in brightness. Indeed, the outburst of BL Lac 
2000 appears as a modest event when compared to the 
1997, 2001, and 2002 outbursts. It is evident that these 
main outbursts show a chromatic behaviour, even if the 
colour changes observed during the short-term flares are 
relatively stronger. By looking at the 2001 outburst one 
can see that the total range covered by the B — R index is 
about twice the range covered during fast flares, whereas 
the ratio between the corresponding magnitude variations 
is about three. This situation is better displayed in Fig. El 

The long-term achromatism found by Villata et al. 
( 2002) in BL Lac 2000 appears to be due to a bluer-than- 
usual pre-outburst state, which we can now recognize as 
a peculiar event. Hence, from now on we will speak of 
"strongly chromatic" fast flares and "mildly chromatic" 
long-term variations. 

We now investigate whether the global behaviour of 
the source can be uniformly described in terms of these 



two components. We try to decouple the two components 
in a way similar to that followed by Villata et al. ( 2002 ) . 

In the upper panel of Fig. we have plotted the B — R 
colour indices versus the galaxy- removed R magnitudes. 
Red (grey) dots represent the data of BL Lac 2000, litera- 
ture data are displayed as light cyan (light grey) dots, and 
data from the BL Lac 2001 campaign plus extension are 
plotted with different symbols according to the observing 
period (see caption for details). One can see the different 
behaviour during the two campaigns: while B — R data of 
BL Lac 2000 are located along two separated quasi-linear 
trends with similar slopes and colour ranges (thus reflect- 
ing the quasi-achromatic long-term behaviour), the other 
ones are more uniformly distributed. As a whole, the in- 
dex dataset is mostly comprised within a parallelogram 
with slopes 0.1 and 0.4. 

A 0.1 slope trend means a mildly chromatic behaviour 
that we have recognized to be typical of the long-term vari- 
ations. On the contrary, a 0.4 slope trend implies strongly 
chromatic variations, as found for the fast flares. The BL 
Lac 2000 points near the upper-left corner of the paral- 
lelogram, which represent data from the pre-outburst, are 
essentially aligned with the steeper side, but their colour 
appears to be bluer than average. Thus, apart from this 
discrepancy, we are presenting a scenario where fast flares 
imply moving along the steeper direction in the B — R 
versus R plot, the position inside the parallelogram de- 
pending on the source brightness state. In other words: 

B-R= constant + 0A(R - R spi ) + 0.1R spi , 

where i? sp i is the brightness level of the long-term compo- 
nent. 
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Fig. 7. B — R colour indices versus R magnitudes after the host-galaxy removal (upper panel): BL Lac 2000 and 
literature data are represented by red (grey) and light cyan (light grey) dots, respectively, while BL Lac 2001 extended 
data are plotted with green (grey) circles (BL Lac 2001 extension), blue (dark) squares (best-sampled period of BL 
Lac 2001, the same of Fig.|B}, and cyan (grey) crosses (remaining part of BL Lac 2001); most of the data lie within a 
parallelogram whose sides have slopes of 0.1 and 0.4. In the bottom panel the same data have been corrected for the 
long-term trend as explained in the text. 



The steeper 0.4 slope has been inferred from the very 
well-sampled data of the last week of the core campaign 
2000, where we can assume that the long-term component 
does not change significantly. Indeed, by looking at Figs. 
7 and 8 of Villata et al. (20|2J), one can estimate d(B — 
R)/dR w 0.4, and a linear correlation analysis gives a 
slope of 0.383 ± 0.004. 

The value of the 0.1 slope has been determined as il- 
lustrated in the next subsection. 



4.1. Removal of the long-term trend 

As done in BL Lac 2000, we performed a cubic spline in- 
terpolation across the binned i?-band flux light curve. To 
have a good interpolation through very fast and large vari- 
ations, like those occurring during the 2001 and 2002 out- 
bursts, we adopted a 4-day binning. The resulting spline is 
plotted as a green (grey) line in Fig.|S| showing the folded, 
galaxy-subtracted flux light curve. 

To remove the long-term trend from the light curve, 
we have rescaled each original flux by dividing it by the 
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Fig. 8. Galaxy-subtracted i?-band flux light curve from 1997 to 2002 (folded, symbols as in Fig.Q; the green (grey) 
line represents a cubic spline interpolation through the 4-day binned light curve; the blue (dark) line, obtained by 
normalizing the spline to its minimum value, represents the actual correction factor for removing the long-term trend. 



figure9.jpg 



Fig. 9. B — R and i?-mag curves as obtained after correcting for the long-term trend (symbols as in Fig.^ horizontal 
lines as in Fig.|SJ); the green (grey) line at R = 16.061 in the bottom panel represents the minimum level of the spline 
(see Fig. |HJ to which the light curve has been normalized. 



ratio between the value of the spline at the considered time 
and its minimum value, C^it) = [F sp i(t)/F m i n )ji (where R 
stands for R band). This correcting factor is shown in Fig. 
|H1 as a blue (dark) line highlighted by the yellow (shaded) 
area. 

Because of the assumption of an achromatic long-term 
trend, in BL Lac 2000 the £?-band flux light curve was 
rescaled by the same time-dependent factor. Now, the cor- 
recting factor for the B fluxes must be found. Instead of 
tracing a cubic spline interpolation across the B fluxes, 
which would be affected by a different and worse sam- 
pling, we search for the best-fit relationship between the 
long-term trends in the two bands. Then we suppose that 
Csit) is some power of Cn(t). Various best-fit procedures 
on different data sub-samples have been applied to find 
the value of the power index; this value has been found 
to range between 1.08 and 1.12. We adopted 1.1 which, 
passing from fluxes to magnitudes, explains the 0.1 slope 
of the B — R versus R plot. 

Once both the R and B fluxes have been corrected 
for their long-term trends, the resulting colour indices are 
shown in the bottom panel of Fig. What we expect to 
see in this plot is the signature of the fast flares alone, 
shifted to the minimum brightness level of the long-term 
trend (R — 16.061). The final result is fairly satisfactory, 
but one can still see some inhomogeneous behaviour, like 
the expected bluer level of the BL Lac 2000 pre-outburst 
data. A further insight into this result is given by Fig. 



El where the residual R magnitudes (bottom panel) and 
corresponding B — R colour indices (upper panel) have 
been plotted versus time. By comparing this figure with 
Fig. one can see that the outburst signatures have been 
completely smoothed away in the i?-band light curve, as 
expected, and almost completely smoothed away in the 
B — R plot, with the main exception of the bluer-than- 
average BL Lac 2000 period. Moreover, one can also notice 
that the variability range of the fast flares in the bottom 
panel is fairly constant around 1 mag, but it seems to 
exceed 1.5 mag or more during the 1997 outburst. This 
represents another peculiarity of an outburst with respect 
to the others. Part of the amplitude excess may be ex- 
plainable by inhomogeneity of the literature datasets, but 
not all of it. 

5. Conclusions 

We have presented U BVRI light curves of BL Lacertae 
from 1994 to 2002 composed from data taken by the mem- 
bers of the WEBT collaboration as well as literature data. 

The colour analysis we performed on the best-sampled 
period of these light curves (1997-2002) leads to the fol- 
lowing results. 

The variability observed in the optical light curves can 
be interpreted in terms of two components: a "mildly- 
chromatic" longer-term component and a "strongly- 
chromatic" shorter-term one. The longer-term component 
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has a typical time scale of a few days, while the shorter- 
term one is responsible for the very fast, intraday flares, 
and determines a steep bluer-when-brighter slope of about 
0.4 in the B — R versus R plot. Moreover, our analysis 
suggests that, in the long-term, the S-band flux level is 
related to the i?-band one according to a power law with 
index 1.1. This means that long-term variations trace a 
0.1 slope in the B — R versus R plot. In other words, all 
possible flux variations fill a parallelogram with 0.1 and 
0.4 slopes in the B — R versus R plane. 

The analysis presented in this paper shows that, in 
general, the BL Lac optical behaviour can be uniformly 
described in terms of these two components fairly well, 
with the caution that the presented decoupling method 
can work well only during well-sampled periods, where the 
light-curve binning can actually represent the long-term 
trend. However, peculiar events also occurred, as in the BL 
Lac 2000 period, and our scenario seems to be more a good 
compromise between slightly different behaviours rather 
than a definitive solution. Other factors appear to change 
in the source history, in one or both of the components. 

Changes in the flaring component could explain the 
larger amplitudes of the 1997 outburst and the bluer BL 
Lac 2000 period, but the remaining blue residues of the 
1997, 2001, and 2002 outbursts suggest other possibilities. 

In Villata et al. (2002) the quasi-achromatism of the 
long-term variations was interpreted in terms of Doppler 
factor variations. If the intrinsic source spectrum around 
the infrared-optical bands is well described by a power law 
(F v oc v~ a ), a Doppler factor variation does not imply a 
colour change. What about a mildly-chromatic behaviour? 
It could be due to a Doppler factor variation on a spectrum 
slightly deviating from a power law (convex for a bluer- 
when-brighter trend). If the spectrum shape changes in 
time, the chromatism of the outbursts also varies, and the 
0.1 B — R slope adopted in this paper might represent only 
a mean value among different slopes. In other words, any 
period should be treated with a proper slope. A further 
investigation of this is beyond the scope of this paper. 

The strongly-chromatic fast flares are likely due to 
intrinsic phenomena, such as particle acceleration from 
shock-in-jet events, widely described in the literature (see 
e.g. Mastichiadis & Kirk 2002 lor a review). Other mod- 
els are also discussed in the recent paper by Bottcher & 
Reimer l|2TM|l . 

Thus, the present extension/enrichment of the BL Lac 
optical dataset partially confirms previous results leading 
to a more refined scenario, but also gives rise to new in- 
triguing questions. 
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